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Upon the introduction of a gas bubble into a liquid possessing a uniform 
thermal gradient, an unsteady thermocapillary flow begins. Ultimately, the 
bubble attains a constant velocity. This theoretical analysis focuses upon 
the transient period for a bubble in a microgravity environment and is 
restricted to situations wherein the flow is sufficiently slow such that 
inertial terms in the Navier-Stokes equation and convective terms in the 
energy equation may be safely neglected (i.e., both Reynolds and Marangoni 
numbers are small). The resulting linear equations were solved analytically 
in the Laplace domain with the Prandtl number of the liquid as a parameter; 
inversion was accomplished numerically using a standard IMSL routine. In the 
asymptotic long-time limit, our theory agrees with the steady-state theory of 
Young, Goldstein, and Block. The theory predicts that more than 90 percent of 
the terminal steady velocity is achieved when the smallest dimensionless time, 
i.e., the one based upon the largest time scale — viscous or thermal — equals 
uni ty. 
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1 . INTRODUCTION 


Forces other than gravity are expected to often dominate the migration of 
bubbles in outer space. In particular, thermocapillary forces will cause 
bubbles to migrate In a thermal gradient. Consider the bubble depicted In 
Fig. 1. The Interface next to the hot liquid will be hotter than that next to 
the cold liquid. Assuming that surface tension varies Inversely with 
temperature, as Is often the case, the cold Interface will have a higher 
surface tension than the hot Interface. The thermocapl 1 lary stresses are such 
that near the Interface liquid flows from the hot side to the cold side, as 
indicated by the arrows. Viscous stresses within the fluid oppose these 
thermocapillary stresses, causing the bubble to migrate toward the hot side. 

If the thermal gradient Is constant throughout the region, the analysis of 
Young et al . (1959) revealed that the bubble will ultimately migrate with a 
steady velocity. 

In the processing of materials In outer space, nonlsothermal bubble-liquid 
systems will be common. Control of these systems will require a good 
understanding of thermocapillary migration. For example, bubbles form In 
glass melts during the manufacturing process; their removal Is essential for 
the glass to be useful. On Earth, gravitational forces aid In their removal. 
During the containerless processing of glass In the microgravity environment 
of space, thermocapl 11 ary- forces are expected to assist their removal (Mattox 
et al . 1982, Subramanlan 1981). Thermocapillary migration may also be 
Important In the design of two-phase heat exchangers for use In outer space. 

A poor design may unexpectedly permit bubbles to migrate en masse to the 
heating surface. A layer of gas next to the hot surface would act as a thermal 
Insulator and prevent the efficient transfer of heat to the liquid phase. 

In contrast to prior experimental and theoretical studies on 
thermocapillary migration (Balasubramanlam and Chal 1987, Hardy 1979, Mattox 
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et al. 1982, Subramanlan 1981, Szymczyk et al . 1987, Young et al . 1959), which 
focus upon steady systems, we here consider the unsteady development of 
thermocapillary migration. Analysis of unsteady thermocapl 1 lary migration Is 
Important for several reasons. First, to study the steady state, one must be 
able to estimate the time required for transients to die out. Also, some 
unexplained results of supposedly steady-state experiments (Neuhaus and 
Feuerbacher 1986, Slekmann et al . 1986), e.g., velocities less than predicted 
steady-state values, may possibly be explained by an unsteady analysis. 

Lastly, the transient period Is expected to be long for relatively large 
bubbles; a theory for their transport will be useful. 

2 . FORMULATION 

Consider a thermocapl 1 lary experiment conducted In outer space as shown 
In Fig. 2. An enclosure consisting of end plates maintained at different 
temperatures and of Insulated side walls Is filled with a liquid. In the 
absence of buoyancy effects, there Is no fluid motion and, as depicted In the 
figure, the temperature field becomes linear when thermal equilibrium Is 
obtained. The experiment then begins with the Introduction of a small gas 
bubble. We neglect the small Immediate effect of the bubble's Introduction 
upon the velocity and temperature fields. Soon thereafter the conduction of 
heat to the Interface causes thermocapillary stresses to develop and migration 
to begin. As Is commonly done, we assume that both viscosity and density of 
the liquid phase are essentially constant throughout the volume, that 
transport processes within the bubble phase may be neglected relative to those 
In the liquid phase, and that surface tension Is a linear function of 
Interfacial temperature. To be definite, we assume that oj - do/dT Is 
negative such that migration is toward the hotter fluid. We also assume that 
the bubble remains spherical during the acceleration; this Is subject to a_ 
posteriori verification. Finally, we assume that both fluid Inertia and the 
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convective transport of heat may be neglected. These latter assumptions 
restrict the analysis to low Reynolds and Marangonl (thermal Peclet) numbers. 

Given the above assumptions, the dimensionless temperature field T 
obeys the unsteady energy equation 

p r = 7 2 T (1) 

In the region outside of the bubble. Here, the liquid's Prandtl number 
Pr ■ v/a, the ratio of kinematic viscosity to thermal dlffuslvlty, appears 
multiplying the left-hand side. The variable t Is time divided by the 
viscous time scale a 2 /v, with the characteristic length scale "a" being the 
radius of the bubble. The dependent variable T Is the scaled difference 
between the actual temperature and that undisturbed temperature which prevails 
In the plane of the bubble's center. The characteristic temperature scale Is 
the product aA, where A Is the undisturbed temperature gradient In the 
system. Eq. (1) Is written with the bubble’s center as the origin of the 
coordinate system. 

Similarly, the velocity field obeys the unsteady Stokes equation 
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In a reference frame moving with the bubble. The vector v Is the scaled 
difference between the fluid's and bubble's velocities. The velocity scale 
for both v and U « kll(t), the scaled bubble velocity. Is the positive 
quantity -OyaA/y. Because the bubble's reference frame Is nonlnertlal, there 
appears In Eq. (2) the fictitious force dU(t)/dt. This force adds to the 
gradient of the hydrodynamic pressure p and vanishes when a steady bubble 
velocity Is obtained. 

At the Interface, viscous and thermocapl 1 lary stresses balance at each 
point. Now both T and v are axl symmetric with respect to an axis that 
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passes through the bubble's origin and Is parallel to the unit vector k, 
which points In the direction of the temperature gradient. It follows that 
T * T(r,e,t> and v = v(r,e,t> and that the stress balance is given by 


8 _ 
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( 3 ) 


Here, r is the dimensionless radial coordinate and 0 Is the polar angle; 
these are depicted in Fig. 2. Since the flow Is axlsymmetrlc, the vector v 
has only two nonzero components: v r and vq, the velocities In the radial 

and polar directions, respectively. 

Because we neglect bubble Inertia and gravitational forces, the total 
hydrodynamic force F h on the bubble vanishes for all t > 0: 

2ir IV 

d<t> n • P <l,0,t)s1n 0 d0 - 0 (4) 


F h (t) - 


JO 


JO 


Here, n Is a unit vector normal to the spherical bubble surface and P<r,0,t) 
Is the hydrodynamic pressure dyadic (tensor). 

We briefly mention the remaining conditions on the temperature and 
velocity fields. The undisturbed temperature and velocity fields prevail (1) 
Initially for all space outside the bubble and (2) far from the bubble for all 
time. The neglect of thermal transport within the bubble requires a zero 
thermal flux normal to the bubble surface for t > 0. Because the bubble does 
not change In size, the radial velocity v r vanishes at the Interface. 

Lastly, the assumption of Incompressibility requires V • v « 0. 


3. SOLUTION 


To solve this system of equations, we defined a modified pressure field 
such that Its gradient is given by the quantity within parentheses In Eq. (2), 
Introduced the streamfunctlon for the axlsymmetrlc flow, and then applied the 
Laplace transform technique to replace the variable t with the parameter s. 
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Analytic functions were found for the transformed temperature and 
streamfunctlon fields that satisfy the differential equations and all 
conditions. The solution yielded the following expression for the transformed, 
dimensionless bubble velocity 0(s): 


t ' * 


U(s) 


1 + VsPr \ 
+ 2 VsPr + s Pry 


+ s(3 fVT) 


(5) 


6(1 +Vs>J 

Multiplication of this expression by s and taking the limit as s 
approaches zero yields a dimensionless terminal velocity of 1/2, In agreement 
with the analysis of Young et al. (1959). We also applied the Initial-value 
theorem twice to Eq. (5) to determine dU/dt at t - 0. The first use of the 
theorem gave the expected result that U(0) * 0. The second use gave 

6 . ( 6 ) 


dU 

dt 


t«0 


The step-function character of the acceleration precluded further application 
of the Initial-value theorem. 

To confirm that the bubble remains spherical throughout the acceleration 
period, we applied our analytic result for the transformed stream function to 
the transformed normal stress condition. The exact satisfaction of this 
condition Indicates that there Is no deviation from the spherical shape. 

Numerical Inversion via an IMSL routine gave the bubble velocity as 
functions of time with Prandtl number as a parameter. Figure 3 Is a log-log 
plot of the dimensionless velocity versus dimensionless time for Prandtl 
numbers of 0.01, 1.0, and 100. 


4. DISCUSSION 

The two most Interesting features of Fig. 3 are that (1) all three curves 
are nearly linear before a transition to the terminal velocity and (2) the 
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dimensionless time to the terminal velocity ranges from order one to order 
100, depending upon the Prandtl number. These features are discussed below. 

For t << 1 , the dimensionless velocity IKt) is given 
approximately by 

U(t) ~ 6t. (7) 

The value of 6 (cf Eq. (6)) derives from the acceleration calculated from 
application of the initial value theorem. Equation (7) is represented in 
Fig. 3 as the dashed line; it fits the data reasonably well for small values 
of t. 

Written in dimensional form, Eq. (7) becomes 



for t*v/a^ << l. (The raised asterisks denote dimensional quantities.) 
According to Eq. (8), the leading term in the short-time expansion for the 
dimensional velocity is independent of viscosity. This result is analogous to 
that for a drop or solid sphere accelerating from rest due to buoyancy forces 
(Chisnell 1987). Though viscosity does not appear explicitly in Eq. (8), it 
does appear in the time scale and thus affects the values of t* for which 
the relation is valid. For highly viscous liquids or small bubbles, the 
relation is valid only briefly. On the other hand, for low-viscosity liquids 
or large bubbles, Eq. (8) is likely to apply for a much longer period of 
time . 

Another interesting feature of Eq. (8) is that the bubble's acceleration 
is inversely proportional to its radius. Consider two bubbles of different 
radii released simultaneously into the same system. During the time regime 
for which Eq. (8) is applicable for both bubbles, the smaller bubble may be 
expected to move faster than the larger one, even though the terminal velocity 
of the smal ler bubble. 
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-o^aA 
* ~ ’ 

Is less. 

Our focus now shifts to estimating the duration of unsteady migration. 
According to Fig. 3, the time t to a terminal velocity Is of order unity If 
the Pr < 1 , but of order 100 If Pr = 100. (Recall that the factor used to 
nondlmenslonal Ize time was the viscous scale a 2 /v.) It thus appears that for 
Pr < 1, the viscous scale Is appropriate, but for Pr > 1 , the product of the 
Prandtl number and the viscous scale Is appropriate. This product equals 
a 2 /a, which Is the characteristic thermal time scale for the system. Observe 
that for any value of Pr, the larger of the viscous or thermal scales appears 
to be a good estimate of the transient period. 

The above rule to predict the time to steady state was based upon 
graphical data for three curves. For a more quantitative picture, we examined 
numeric data for five liquids having Prandtl numbers In the range 0.01 < Pr < 
100. For each Prandtl number, we determined the dimensionless times t go and 
x g5 (based upon the appropriate scale) for a bubble to attain 90 and 95 percent 
of the terminal velocity. All results (Table 1) are of order unity and thus 
confirm that the rule Is correct. 

To understand why the above scaling works, observe that If Pr << 1, the 
thermal field attains Its steady-state value while the fluid motion Is just 
beginning to develop. The time for the velocity field to reach Its steady 
state Is thus governed by the viscous time scale. If Pr >> 1, the fluid 
motion, which Is driven by thermocapl 1 lary forces and thus by the temperature 
field, becomes quasistatic with a slowly-developing thermal field. In this 
case, It Is the thermal time scale that determines when the bubble reaches Its 
terminal velocity. 
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5. SUMMARY 


The unsteady thermocapl 1 1 ary migration of a bubble within an otherwise 
quiescent liquid was studied via analytical methods. Within the same system, 
small bubbles are predicted to Initially move more rapidly than larger ones, 
even though their terminal velocities are less. The larger of the viscous and 
thermal time scales Is a good estimate of the duration of the unsteady 
migration period. 
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TABLE 1. - DIMENSION- 
LESS TIMES TO TERMINAL 
VELOCITY 


Pr - v/a 

t 90 

x 95 

0.01 

0.28 

0.47 

.10 

.30 

.52 

1.0 

.58 

1.10 

10.0 

.29 

.68 

100.0 

.26 

.65 




IFO T = ^f < 0, ... 




FI6URE 1. - ORIGINS OF THERMOCAPILLARY 
MIGRATION. VARIATIONS IN SURFACE TEN- 
SION ALONG A NONISOTHERMAL INTERFACE 
CAUSE THERMOCAPILLARY SHEARING STRESSES, 
WHICH VISCOUS SHEARING STRESSES TRANSMIT 
TO THE LIQUID. LIQUID IMPOSES A RE- 
ACTION FORCE UPON THE BUBBLE, CAUSING 
IT TO SWIM TOWARDS THE HOT LIQUID WITH 
A VELOCITY U(t). 


FIGURE 2. - SCHEMATIC OF A THERMOCAPILLARY EXPERIMENTAL CELL IN ZERO 
GRAVITY. IN TYPICAL EXPERIMENTS, A LINEAR TEMPERATURE FIELD IS 
ESTABLISHED PRIOR TO THE INTRODUCTION OF A GAS BUBBLE AT POINT 0. 
THERMOCAPILLARY ACTION CAUSES THE BUBBLE TO MIGRATE WITH THE UN- 
STEADY VELOCITY U(t) IN THE DIRECTION OF THE THERMAL GRADIENT. 
EVENTUALLY, THE MIGRATION VELOCITY BECOMES STEADY. 



FIGURE 3. - DEVELOPMENT OF THERMOCAPILLARY MIGRATION. 

EACH SOLID CURVE REPRESENTS A PLOT OF THERMOCAPILLARY 
MIGRATION VELOCITY U (t) VERSUS TIME t FOR A SPECI- 
FIED PRANDTL NUMBER PR = V/a. THE DASHED LINE PROVIDES 
A GOOD ESTIMATE OF THE VELOCITY FOR SMALL VALUES OF t. 
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